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Evidence for neutrino oscillations is pointing to the existence of tiny but finite neutrino masses. Such 
masses may be naturally generated via radiative corrections in models such as the Zee model where 
a singlet Zee-scalar plays a key role. We minimally extend the Zee model by including a right- 
handed singlet neutrino vr. The radiative Zee-mechanism can be protected by a simple U(l)x 
symmetry involving only the vr and a Zee-scalar. We further construct a class of models with a 
single horizontal (7(1)fn (a la Froggatt-Nielsen) such that the mass patterns of the neutrinos and 
leptons are naturally explained. We then analyze the muon anomalous magnetic moment (g^ — 2) 
and the flavor changing /1 — > ey decay. The ur interaction in our minimal extension is found to 

induce the BNL g M — 2 anomaly, with a light charged Zee-scalar of mass 100 — 3 00 GeV. 
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The standard model (SM) of the electroweak and 
strong interactions has to be extended, in light of the ex- 
isting neutrino oscillation data jjj which provides strong 
evidence for tiny but finite neutrino masses. The squared 
neutrino-mass differences are found to have two distinct 
ranges, 1CT 11 eV 2 < 8m% < lCT 5 eV 2 and <5m 2 tm ~ 
10 -3 eV 2 , in order to interpret the solar and atmospheric 
neutrino anomalies, respectively. The very small neu- 
trino masses may be explained by invoking the seesaw 
mechanism [|| where a heavy right-handed Majorana sin- 
glet neutrino is introduced at the grand unification scale 
(10 10-16 GeV). On the other hand, the radiative mech- 
anism for neutrino mass generation, such as advocated 
in the Zee model ||, provides an important alternative, 
where the relevant new physics is expected to show up at 
or near the weak scale. Such low scale models, besides 
being able to explain some features of neutrino oscilla- 
tions such as the bi-maximal mixing [Q, are clearly of 
phenomenological relevance. The existing precision data 
have put some nontrivial constraints on the models 
and further tests will be available at forthcoming collider 
experiments. 

The minimal Zee model || contains the three active 
left-handed neutrinos of the SM and a bilepton singlet 
Zee-scalar which plays a key role for radiative genera- 
tion of their Majorana masses. There is no underlying 
reason that forbids the existence of light right-handed 
singlet neutrinos (vr), as vr can also be naturally con- 
tained in various extensions of the SM (such as the mod- 
els with left-right symmetry or 50(10)). The introduc- 
tion of vr thus provides the simplest possible extension 
of Zee model. However, a simple embedding of vr results 
in the loss of the prediction of neutrino masses as they 
are rendered arbitrary by the tree-level Dirac mass terms 
iQiPPfll' ^ n * ms wor k! we first build a class of minimally 
extended Zee models including vr (called Type-I) and in- 
voke a simple U(l)x symmetry (or its discrete subgroup 
Z n ) to effectively protect the radiative Zee- mechanism 



by forbidding the mixings between the vr and the ac- 
tive neutrinos. Such an extension is nontrivial since the 
successful embedding of a singlet vr requires the addi- 
tion of a second singlet Zee-scalar in the minimal exten- 
sion. Next we note that the original Zee model neither 
predicts the size of the Zee-scalar Yukawa couplings nor 
provides any insight on generating the lepton masses and 
their hierarchy. Though our simplest U(l)x m Type- 
I models protects the radiative Zee-mechanism, we can 
use this same U(l)x group as a horizontal symmetry 
involving both neutrinos and leptons (called Type-II), 
and thus explain the mass patterns of the neutrinos and 
leptons in a natural way, a la Froggatt-Nielsen |J. In 
both Type-I and -II models, the vr can interact with the 
right-handed muon and Zee-scalar with a natural O(l) 
Yukawa coupling, which has striking phenomenological 
consequences. 

Finally, we apply the Type-I and -II models to ana- 
lyze the Zee-scalar-induced contributions to the muon 
anomalous magnetic moment g^ — 2 and the lepton- 
flavor- violation decay /j, — > e-/ . We find that the re- 
cent BNL g M — 2 anomaly JlO| can be explained with 
a light charged Zee-scalar of mass around 100 — 300 GeV. 
Our models also have the fi — > decay branching ratio 
around or below the current experimental limit. 

Minimal Extension of the Zee Model 
with a Right-handed Singlet Neutrino 

The minimal Zee model || introduces one extra singlet 
charged scalar (Sf) together with the usual two-Higgs- 
doublet sector. By assuming no right-handed vr, as in 
the SM, this scalar only interacts with left-handed neu- 
trinos and leptons. Thus, the Zee model contains the 
following additional Lagrangian, 
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/as (^T L -^ L ) + m 3 (tff ^-^i/")] S++h.c. 

where € (e, /i, r) and Lj = (i/j, Ijl) T is the left- 

handed doublet of the jth family. The (f/i, H 2 ) are the 
usual two-Higgs-doublets with hypercharge (1/2, —1/2), 
where H x = {-H+ , H?) 71 , ff 2 = (#2- ff 2 ~) T , and 
Hj = ir 2 H* . The Yukawa sector can conserve total lep- 
ton number by assigning to Sf 1 the lepton numbers +2. 
Thus, the total lepton number is only softly violated by 
the dimension-3 trilinear Higgs operator in Eq. (Q). As 
such, the small Majorana neutrino masses are radiatively 
generated at one-loop and are automatically finite. 

We minimally extend the Zee model by including a 
single right-handed Dirac neutrino vr with following 
Yukawa interactions, 

AC 2 = [fi^e R + f 2 ~i^tiR + f3~t R ~T R \S+ + h.c. (2) 

where S^ is a second singlet Zee-scalar. The nontrivial 
issue with embedding vr is to avoid arbitrary tree-level 
Dirac mass terms generated by the Yukawa interactions 
LjH 2 vr and LjH\Vr (which mix Vj and vr), so that 
the predictive power of the radiative Zee-mechanism can 
be effectively protected. We achieve this goal by not- 
ing that the Yukawa sector (j|) of vr possesses a global 
U(l)x symmetry, which can properly forbid the neutrino 
Dirac mass terms once a Zee-scalar is included to- 
gether with vr. It can be shown, by assigning the most 
general U(l)x quantum numbers for the Zee-model with 
Vr, that the unwanted tree- level neutrino Dirac masses 
cannot be removed without ■ We define our simplest 
Type-I models with U(X)x in Table 1, where only vr and 
carry XJ{V)x charges while all other fields are singlets 
of U{V)x- Hence, the Type-I extension gives a truely 
minimal embedding of vr into the Zee model. 



Table 1. Quantum number assignments for Type-I and -II 
models. The hypercharge is denned as Y = Q — I3 . 
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Table 2 classifies all (dis-)allowed operators of Type- 
I up to dimension-4. It shows that, as long as x ^ 0, 
the radiative Zee-mechanism is protected and the vr 
remains massless. Such a massless vr does not con- 
tribute to the invisible Z-width as it carries no weak 
charge. A special case of our Type-I is to consider its 
discrete subgroup Z4 under which vr and S~i transform 
as, vr — » wr , S2 — > TiS^ 1 , while all other fields remain 



invariant. Other non- minimal variations of our Type-I 
can be easily constructed. 



Table 2. Summary of U(l) charges carried by the effective 
operators in Type-I and -II models. 
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Neutrino Oscillations, Lepton Masses and 

Horizontal [/(1)fn Symmetry 

While the above Type-I models give the most economic 
embedding of vr with all the good features of the original 
Zee-model retained, they do not provide any insight on 
two important issues: (i) There is no theory prediction on 
the size of the Zee-scalar Yukawa couplings fjf in Eq. (I), 
but the neutrino oscillation data requires the following 
hierarchy pj: 

^4-3xl0 2 , £» * * 10 2 or 10 7 , 

Ji3 f 23 6m z Q ^ 

where /13//23 — 10 2 (10 7 ) corresponds to the MSW large 
angle solution (vacuum oscillation solution), (ii) The 
small lepton masses and their large hierarchy are not un- 
derstood. Our goal is to construct this same U(l) group 
as a horizontal symmetry involving all the leptons so that 
these two issues can be naturally explained a la Froggatt- 
Nielsen (FN) @. [This U(l) will be called U(1) FN .] 
The basic idea is to consider a horizontal J7(1)fn sponta- 
neously broken by the vacuum expectation value (5°) of 
a singlet scalar S°. We can assign ?7(1) FN charges for rel- 
evant fields such that different mass terms are suppressed 
by different powers of e = (S°)/A where A is the scale 
at which the £7(1)fn breaking is mediated to the light 
fermions. For instance, a low energy effective operator 
carrying a net £7(1)fn charge q (either > or < 0) will 
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be suppressed by e' 9 '. Though all mass terms are now 
allowed in the effective theory, we will build a class of 
FN-type models (called Type-II) in which the arbitrary 
tree-level neutrino Dirac-mass terms are suppressed to 
a level much below the one-loop radiative Zee-masses, 
and thus the predictive power of the Zee-mechanism re- 
mains. The role of the FN-scalar S° is to provide the 
spontaneous Z7(1) FN breaking and generate the relevant 
f^(l)FN-invariant effective operators that will give the de- 
sired neutrino Yukawa couplings and lepton masses at the 
weak scale. The heavy S° will be eventually integrated 
out from the low energy theory and our relevant particle 
spectrum of Type-II is the same as Type-I. 

We provide two typical Type-II constructions, called 
Type-IIa and -lib, respectively. The Type-IIa is the sim- 
plest extension of Type-I by further involving only the 
right-handed weak-singlet leptons in the £/(1)fn (cf. Ta- 
ble 1). In the Type-lib models, we further assign each 
lepton doublet Lj a charge Uj. So, the lepton masses are 
determined by IjR charges in Type-IIa, while Type-lib 
determines these masses by the charges of both £jR and 
Lj . The low energy effective operators up to dimension-4 
(with the heavy S° integrated out) are classified in Ta- 
ble 2, from which we derive the general conditions for 
protecting the Zee-mechanism in Type-II, 

10 > ~ |a;| > 1, and \x - uj\, \x + y\ > 1, (4) 

with xx' > and \y\, \z\ ~ 0(1). For the explicit anal- 
ysis below, we choose a typical value of the suppression 
factor e ~ 0.1. Thus, choosing leptons in mass-eigenbasis, 
we write their mass ratios as, 



m t 



e 4 : e 1 : e° . 



(5) 



which require, 

(yi-wi)-(y 3 -"3)=±4, (2/2—"2)- (2/3-1*3) = ±1- (6) 

The tau Yukawa coupling itself can be estimated as 
y T ~ (m T /mt) tan/3 ~ 10~ 2 tan/3 ~ e 1 (with tan/3 = 
(Hz)/ (H\)), in the typical range of tan/3 ~ 10 — 40, and 
this restricts the C/(1)fn charges of r as y 3 — u 3 = ±1. 
Table 3 summarizes three explicit realizations of Type-II 
models. From Table 3 and Eq. (||), the Yukawa couplings 
of vr are predicted as, 

Type Ila : (f u f 2 , f 3 ) ~ (e 3 , 1, e 1 ) ; 
Typellbl: (/1, / 2 , fa) ~ (e 5 , 1, e 3 ) ; (7) 
TypeIIb2: (f u f 2 , / a ) ~ (e 10 , 1, e 3 ) . 

From Table 3 and Eq. (|l|) , we further predict the left- 
handed Yukawa couplings fjji , 



Type Ila 
Type Ilbl 
Type IIb2 



(/12, fi 3 , /2s) ~ 5 

(/i2,/i3, / 23 )^(e 4+z ,e 6+z ,e 8+2 ); 

(/i2,/i3,/ 23 )^(e 3+2 , e 5+z ,6 12 +-): 



(8) 



where the allowed values of z are defined in Table 3. 
Thus, Type-IIa suppresses fjji couplings to O(10~ 2 — 
10~ 4 ). The models in Type-IIbl (-IIb2), however, nicely 
accommodate the hierarchy (||) for the MSW large an- 
gle solution (vacuum oscillation solution), while the pre- 
dicted size of fi2 ~ 10~ 3 — 10~ 6 is also of the right order 
H . Finally, it is trivial to extend these models with more 
than one singlet Vr (i.e., VRj with j — 1, • • • ,N UR and 
N VR — 3 for instance), by simply defining them to share 
the same U(l) charges as in Tables 1 and 3. 

Table 3. C/(1)fn quantum number assignments for Type-IIa, 
-Ilbl and -IIb2 models. [The defined range of z is \z\ ~ 3 for 
Type-IIa and < 2 < 3 for Type-IIbl and -IIb2.] 





Li 


Li 


L 3 


eR 




TR 


VR 


Hi 


H 2 


Si S2 


Ila 











-5 


-2 


-1 


x+1 





z 


— z 1 — x 


Ilbl 


-1 


-3 


-5 


4 


-1 


-4 


X 





z 


— z 1 — x 


IIb2 


2 


-5 


-7 


7 


-3 


-6 


X 





z 


—z 3 — x 



Zee Scalars, Muon g — 2 and /i — > e'y 

The above minimally extended Zee-type models eco- 
nomically incorporate the vr and naturally explain the 
mass patterns of the neutrinos and leptons. The Zee- 
scalar Yukawa couplings with the neutrinos/leptons also 
exhibit an interesting spectrum. Now we are ready to an- 
alyze their phenomenological impact. The Brookhaven 
E821 collaboration has announced a 2.6 standard devi- 
ation in the muon anomalous magnetic moment = 



(5u-2)/2, i.e., Aa M = a^P- a « M = (42.6± 16.5) x 10 
p0|, which gives a 90%C.L. range for new physics, 



-10 



15.5 x 10 



-10 



< Aa p < 69.7 x 10 



-10 



(9) 



Different authors ]Tl| ] have interpreted this anomaly in 
terms of supersymmetry, muon compositeness, extra Z' , 
leptoquarks and extended neutrino models. We attempt 
to explain it from the contribution of the Zee-scalars and 
the singlet vr in our minimal Type-(I,II) models. 

The Zee-scalars Sf' and S2 in Type-I/II contribute to 
g M — 2 via the Yukawa couplings /12,2s with (/x/,, v e ,r) 
and / 2 with (fiR, Vr), respectively. Thus, we have, 



m% /|/i 2 | 2 + |/ 23 | 2 , I/2 1 2 



96tt 2 



11.8 x 10" 



x I/2 



Ml 

02 



lOOGeVY 
Ms 2 ) 



(10) 



with Mi = (cos 2 0/M 2 , +sin 2 0/A^, ± )- 2 . Here (M s; , 
M H ,±) are the mass-eigenvalues of the two charged 
scalars in Eq. (1) and 4> is their mixing angle. Our mod- 
els forbid or highly suppress the mixing between 



3 



and Sf [cf^Table2 and Eq. (§)]. Note that the /i 2i23 



terms in (|10|) are negligible compared to the / 2 term 
for Type-I and Type-II (with \z\ > 1), cf. Eqs. @-(|). 
The precision bound from — > v^ev^ decay gives [Q, 
fxifMx < 0.18/TeV, which, combined with Eq. (§), also 
renders the /i 2 ,23 terms irrelevant for the g^ — 2 anomaly. 
Hence, the original Zee-model cannot resolve the g^ — 2 
anomaly. From (||) and ([lO]), we deduce, 



41.1 GeV < -j-^f < 87.2 GeV . 
I/2I 



(11) 



Since the LEP2 direct search for charged particles re- 
quires Ms 2 > 100 GeV, |/ 2 | is constrained as, 



1-1 < IM < O(l) : 



(12) 



where the upper bound is imposed by perturbativity. 
Eq. (|l^) is in the predicted range of / 2 for our Type- 
II models [cf. Eq. (@)]. Thus, combining ([□]) and @, 
we conclude that, to fully accommodate the BNL g^ — 2 
data, the Zee-scalar in our Type-I and -II models has 
to be generically light, around 100 — 300 GeV. This leads 
to the possibility of discovering the light charged Zee- 
scalar at the Tevatron Run-2, the LHC or a future high 
energy Linear Collider. Since / 2 3> /1.3 for Type-II, the 
S% — > [^vr decay has a large branching ratio, and dom- 
inates over the er^VR and t^isr modes. Though Type-I 
models do not predict /i^, they allow / 2 = 0(1) while 
/i,3 ^ /2 is forced by the fj, — > ej and r — > fij bounds 
below. We have the partial decay width, 

r[S± - i4v R ] ~ -^M S2 ~ 0(1) GeV, (13) 



16tt 2 



for Ms 2 — 100 — 300 GeV. Hence, S 2 is a very nar- 
row spin-0 resonance. The predicted branching ratio 
BrfS^ — > H%vr] — 1 in Type- (I, II) suggests that S 2 
can be best detected via muon plus missing energy. 

Our models have further implications for the flavor- 
violating rare decay /i — * ej, whose branching ratio is 
bounded by, Br[^ -> ey] < 1.2 x 10" 11 , at the 90%C.L. 
The Type- (I, II) models give the following contributions, 



Br[/i ej] 



384tt 



I/1/2 1 

Ml 2 



I/13/23 1 2 



Ml 



\hf2\ 2 (v/M S2 ) 4 



(14) 



384tt' 

where v = (V2G F )" 1/2 ^ 246 GeV. Thus, we derive, 

I/1/2 1 < 2.3 x lO" 4 (M S2 /100GeV) 2 . (15) 

Combining this with the — 2 bound in (|lll), we find, 
I/1//2I < 1-8 x lO" 4 , => /!<(2-6)xl0- 4 . (16) 

Comparing this with Eq. (g), we see that Type-IIbl has 
fi just below the current bound while Type-IIb2 is well 
below it. On the other hand, the fi coupling in Type-IIa 



lies slightly above the limit by a factor of 2 — 3; given the 
uncertainty of the parameters, it can be easily adjusted 
to stay within the bound. Also, a much weaker bound 
on /3 can be derived from r — » /17 decay, i.e., I/3I < 
0.06- 0.16 - O(0.1) at 90%C.L., for 1.1 < \f 2 \ < 3, 
which is consistent with the Type-II predictions in (Q). 
Finally, if we include N VR (> 2) singlet i/Rj with the same 
Yukawa coupling the upper [lower] bound in Eq s. (|Tl] ) 
and (|l|) [Eq. @] will be relaxed by a factor of ^JN V ^. 

In summary, the Zee model naturally generates small 
neutrino Majorana masses by radiative corrections, but it 
neither predicts the Zee-scalar Yukawa couplings nor pro- 
vides any insight on the lepton mass hierarchy. We have 
constructed a class of minimally extended Zee-models 
with the right-handed neutrino v-r embedded, where a 
C7(l) symmetry protects the radiative neutrino masses 
while generating the lepton mass hierarchy, the hierarchy 
of the Zee-scalar Yukawa couplings required by the neu- 
trino oscillation data, the hierarchy of Zee-scalar Yukawa 
couplings necessary for consistency with the fj, — > e^f 
bound, and the size of the Zee-scalar Yukawa coupling 
needed for the BNL g„ — 2 anomaly. Furthermore, a 
light Zee scalar S 1 ^ is predicted in our models, with a 
mass around 100 — 300 GeV. 
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